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Abstract 

An individual-tree diameter growth model was developed for Cunninghamia lanceolata in Fujian province, southeast China. 
Data were obtained from 72 plantation-grown China-fir trees in 24 single-species plots. Ordinary non-linear least squares 
regression was used to choose the best base model from among 5 theoretical growth equations; selection criteria were the 
smallest absolute mean residual and root mean square error and the largest adjusted coefficient of determination. To 
account for autocorrelation in the repeated-measures data, we developed one-level and nested two-level nonlinear mixed- 
effects (NLME) models, constructed on the selected base model; the NLME models incorporated random effects of the tree 
and plot. The best random-effects combinations for the NLME models were identified by Akaike's information criterion, 
Bayesian information criterion and —2 logarithm likelihood. Heteroscedasticity was reduced with two residual variance 
functions, a power function and an exponential function. The autocorrelation was addressed with three residual 
autocorrelation structures: a first-order autoregressive structure [AR(l)], a combination of first-order autoregressive and 
moving average structures [ARMA(1,1)] and a compound symmetry structure (CS). The one-level (tree) NLME model 
performed best. Independent validation data were used to test the performance of the models and to demonstrate the 
advantage of calibrating the NLME models. 
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Introduction 

China-fir (Cunninghamia lanceolata (Lamb.) Hook) is the most 
commonly grown afforestation species in southeast China because 
of its fast growth and good wood qualities. It is widely used for 
buildings, furniture, bridge construction and many other purposes. 

Growth and yield models are commonly used for forest 
management planning because they can simulate stand develop- 
ment and production under various management alternatives 
[1,2]. Individual-tree diameter growth models are a fundamental 
component of forest growth and yield prediction frameworks [3- 
5] . The models are based on extensive growth data obtained from 
diverse regions and management levels. Individual-tree diameter 
growth can be expressed as a function of tree size, competitive 
effect, stand structure and site quality [6] . A distance-independent 
individual-tree model structure may be flexible enough to predict 
diameter growth in monospecific even-aged stands and in mixed- 
species and multi-aged stands [7]. 

Regression analysis, such as ordinary non-linear least squares 
(ONLS) regression, is the most commonly used statistical method 
in forest modeling [8]. Individual-tree diameter growth models 
have been fitted to growth increment data collected repeatedly 
over time on the same tree [9] . The hierarchical nature of the data 
results in spatial and temporal correlation among observations 



made in the same sampling unit (i.e., plot and tree) [10]. However, 
the stochastic structure is often ignored and independence of 
observations is assumed [11-15]. Furthermore, the data are 
autocorrelated and cannot be considered independent samples of 
the basic tree population [10]. The ONLS regression assumption 
of independent residuals is therefore violated, biasing the estimates 
of the standard error of the parameter estimates [16]. Many recent 
efforts to develop diameter growth models have used nonlinear 
mixed-effects (NLME) models [5,17]. 

NLME models include both fixed effects, which are parameters 
associated with an entire population or with certain repeatable 
levels of experimental factors, and random effects, which are 
associated with individual experimental units drawn at random 
from a population [18]. Random effects account for spatial and 
temporal correlation by defining the covariance structure of the 
model's random components and by using this structure during 
parameter estimation. NLME models provide an efficient 
statistical method for explicitly modeling hierarchical stochastic 
structure. Growth models can be calibrated by predicting random 
components from tree- or plot-level covariates when a new subject 
is available and is not used in the fitting of the model by using the 
empirical best linear unbiased predictors (EBLUPs) [3,4,12,19]. 

Statistical models in which both fixed and random effects enter 
nonlinearly are increasingly common in the biosciences [20] . The 
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Figure 1. Seventy two trees in twenty four sample plots on Jiangle state-owned forest farm in southeast China. 
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models are relevant to many disciplines, including forestry, 
agriculture, ecology, biology, biomedicine and pharmacokinetics 
[21]. They are used to analyze data with complex structures, 
including grouped data, longitudinal data, repeated measures data 
and multivariate multilevel data [22]. One of the most common 
applications is for analysis of nonlinear growth data [23]; these are 
data measured repeatedly over time on the same tree (multiple 
observations obtained from the same sampling unit or subject in 
sequence over time) and are also known as longitudinal data 
[12,24]. 

The main purpose of this study was to develop an individual- 
tree diameter growth model for C. lanceolata (Lamb.) Hook 
growing in Fujian province, southeast China. The data were 
derived from 144 increment cores from 72 trees in 24 sample plots. 
One-level and nested two-level nonlinear mixed modeling 
approaches that included both fixed and random components 
were applied to the hierarchical structure of the data. This 
diminished the level of variance among the sampling units, which 
were included as covariates. In developing the diameter growth 
models, we considered nested two-level models and a single-level 
model. The first level is the plot and the second level is the tree, 
nested within the plot. Our preliminary analysis showed that the 
NLME models with random effects effectively removed the 
heteroscedasticity and autocorrelations in the repeated-measure 
data and therefore could be important tools for sustainable 
management of China-fir species within the study area. The 
predictive ability of the developed model and the applicability of 
the NLME model were demonstrated using separate validation 
data. 

Materials and Methods 

Data 

The data were obtained from 24 single-species plots of 
plantation-grown China-fir on the Jiangle state-owned forest farm 
in southeast China (Figure 1). One-hundred and forty-four 
increment cores were collected from 72 trees; 15 cores missed 
the pith and were excluded from analysis. The increment cores 
were extracted from three mean trees in each plot; the mean trees 
were trees with a diameter at breast height (dbh; 1.3 m above 
ground) approximately equal to the plot mean dbh. Two cores 
were collected perpendicular to each other from each tree at breast 
height. The sample plots were square and varied in size from 400 
to 600 m 2 . All standing live trees (height >1.3 m) on the plots 



were measured for dbh (outside bark) and tree height. Three to 
five dominant trees on each plot were chosen to calculate plot 
dominant height. Using the Lintab tree-ring measurement system 
(Rinntech Company in Germany), the width of each annual 
growth ring (radial increment data) was measured; trees were 
assumed to be round, so the diameters were calculated as twice the 
radius. Growth data from the seed orchard at the forest farm 
indicate that China-fir requires two years to attain a height of 
1.3 m. Diameter data were therefore assigned an initial age of 3 
years. Independent data were used for model validation. The data 
were randomly divided into two groups; 75% of the points were 
used for model fitting, and 25% were used for model validation. 
The fitting data and the validation data included 54 trees from 23 
plots and 18 trees from 13 plots, respectively. Summary statistics 
for both fitting and validation data are shown in Table 1 . 

Methods 

Nested two-level NLME model. The model data were 
derived from the measured annual increment of the sampled trees. 
The nested sampling structure created a high degree of correlation 
among observations taken from the same tree and plot. The 
mixed-effects modeling approach is a common means of 
addressing the correlation structure in the data [13,23]. A general 
expression for a NLME model can be defined as [22,25]. 

DBH, y/t =f((j> v ,t tik )+e i j k , i=l,...,M,j= U-M, k= l,...,«, y (la) 

where M is the number of plots, M, is the number of trees within 
the ith plot, and fly is the number of observations (increments). 
DBH y 4 is the dbh (cm) at the kth age of thejth tree taken from the 
ith plot, fe* is the age, <j>g is the parameter vector rxl (where r is 
the number of parameters in the model), / is a nonlinear function 
of the predictor variables and the parameter vector, and e,^ is the 
within-group error including the within-group variance and 
correlation [26]; the error is assumed normally distributed with 
a mean of zero and a positive-definite variance-covariance 
structure Ry, generally expressed as a function of the parameter 
vector S [27]. 

s#~N(0,fy) (lb) 
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Moreover, fcj can be expressed as: 
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where X is the p x 1 vector of fixed population parameters (where jf) 
is the number of fixed parameters in the model), fij and l-iy are the 
q\ x 1 and cj2 x 1 vectors of random effects associated with the first 
and second levels, respectively (where q\ and are the numbers 
of random parameters of two-level in the model), which are 
assumed to be normal (or Gaussian) with a mean of 0 and have the 
variance-covariance matrices i//, and ifc; these are the q\ xq\ and 
<?2 x <?2 variance-covariance matrices associated with the first and 
second level random effects, respectively. Ay, _B, ; and By are the 
design matrices rxp, rxq\ and rxq2 for the fixed and random 
effects specific to each sampling unit. 

Individual-tree diameter growth equation. Five theoret- 
ical nonlinear growth equations, the Richards, Weibull, Korf, 
Logistic and Schumacher equations, were selected as candidates 
for modeling diameter growth. These equations are widely used 
for the simulation of individual-tree growth, particularly the 
Richards and Korf equations. Mathematical expressions of the 
equations are shown in Table 2. 

The five above-mentioned equations are all S-shaped growth 
equations with inflection points and asymptotes. A characteristic of 
the Richards, Weibull and Korf equations is that the coordinates 
of the inflection points are variable multiples of asymptotic values; 
in contrast, the equivalent values of the logistic and Schumacher 
equations are fixed multiples [28]. The five equations were initially 
fit by ONLS regression using the R nls function without random 
parameters. Different initial values for the parameters were tried to 
ensure that a global minimum was achieved. The best performing 
function was selected as the base model by applying three 
statistical criteria; absolute mean residual (AMR), root mean 
square error (RMSE), and adjusted coefficient of determination 
(■^adj) P 9 1- The function with the smallest AMR and RMSE and 
the largest R 2 ii: provides the best fit. The adjusted coefficient of 
determination is used similarly as an unbiased estimator in both 
multiple regression and canonical redundancy analysis. The 
formulas of the fit statistics are: 
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where Jyi is the predicted increment at the &th age within the jth 
tree within the zth plot, y is the average of observations. 



PLOS ONE | www.plosone.org 



3 



August 2014 | Volume 9 | Issue 8 | e104012 



NLME Diameter Growth Models 



Table 2. Mathematical expressions of the five equations. 




Equation 


Expression 


Inflection point 




Parameters 






Abscissa 


Ordinate 




Richards 


DBH = 0 1 (l-exp(-0 2 f)) f, 


l/(0 2 ln0 3 ) 


0,(l-l/0 3 )*> 


01. 02>O 


Weibull 


DBH = 0,(l-exp(-0 2 r'-<)) 




#,(I-exp(l-^ 3 )/^ 3 ) 


01, 02, 03>O 


Korf 


DBH=^exp(-^ 2 /(*3) 




^iexp((f,-l)/f,) 


01, 02, 03>O 


Logistic 


DBH=^/(l+catp(fc-^/)) 




0i/2 


01, 03>O 


Schumacher 


DBH = 0,exp(-0,/r) 




0,e~ 2 


01, 02>O 


0i, (j) 2 and 03 are the formal parameters. 
doi:1 0.1 371 /journal.pone.01 0401 2.W02 



Mixed parameter evaluation. A crucial issue in fitting 
mixed-effects models is deciding which parameters should be 
considered random effects and which can be treated as fixed 
effects. A common approach is to start with random effects for all 
parameters and then to examine the fitted object to decide which, 
if any, of the random effects can be eliminated from the model 
[18]. Different combinations of model parameters were therefore 
tested to ascertain their contribution to predictions of diameter 
growth; the best model was selected by Akaike's information 
criterion (AfC) [30], Bayesian information criterion (BIC) [31] and 
— 2 log-likelihood (—2 LL) [32]. The best model gave the smallest 
AIC, BIC and —2 LL. The appropriate variance function and 
autoregressive structure for the NLME models were determined 
by the likelihood ratio test (LRT) [18,33]. All NLME models 
presented in this paper were calibrated using the nine function in 
the R statistical environment [34]. 

Determining the variance-covariance structure. The 
variance-covariance matrices t/'; and lj/y are positive-definite and 
symmetric, which is to say that all their eigenvalues must be strictly 
positive [18]. A hypothetical 2x2 variance-covariance matrix is 
shown as follows [24,35]: 



including both correlation effects and weighting factors. The 
matrix is expressed as [24,36]: 



Rtt 



3 j s-Al- 



(6) 



where for a tree j in plot i, with increment, Ru is the n« xtiy 
within-tree variance-covariance matrix that defines within-group 
variability, Gy is an fly xn v diagonal matrix of within-tree error 
variance (heteroscedasticity), Ijj is an riyXfiij matrix of within-tree 
autocorrelation of the errors, and a 2 is a scaling factor for the error 
dispersion [13]. 

In individual-tree diameter growth models, the variance is often 
found dependent on the means, and the variance will generally 
increase with increasing mean tree diameter. To remove this 
effect, we modeled the variance as an exponential function or 
power function which was used for Gy matrix [18]. And for the 
exponential function and power function, the diagonal elements of 
Gy are & and 28tyk, respectively, and the off-diagonal elements 
are all 0. 



varexp(£ (/ A..) =rrexp(25t ijk ) 



(7) 



2 



where rr^ and g\, are the variance for the random effects u and w, 
respectively, and o m = d wu is the covariance between random 
effects u and w. 

Determining the structure otRy. The matrix Ry is allowed 
to depend on both random and fixed effects, as well as on a set of 
common but unknown parameters. The matrix accounts for 
within-plot heteroscedasticity and autocorrelation [24,26,27] by 



varpower^) 



= ff % 



(8) 



Autocorrelation structures were used for ly matrix to address 
the within-tree autocorrelations of the errors observed in the data 
[37,38]. A method was selected from among three commonly used 
approaches: first-order autoregressive structure [AR(1)], a combi- 
nation of first-order autoregressive and moving average structures 
[ARMA(1,1)], and the compound symmetry structure (CS) [18]. 



Table 3. Performance criteria for individual-tree diameter growth equations. 




Equations 


Fitting data 






Validation data 






AMR 


RMSE 




AMR RMSE 


rU 


Richards 


2.2169 


3.3224 


0.7729 


3.7228 4.7236 


0.7448 


Weilbull 


2.2724 


3.3503 


0.7690 


3.7646 4.7196 


0.7452 


Korf 


2.1286 


3.2710 


0.7855 


3.6101 4.5369 


0.7641 


Logistic 


2.3789 


3.4523 


0.7548 


4.0064 5.0181 


0.7120 


Schumacher 


2.1390 


3.2808 


0.7785 


3.5979 4.6037 


0.7576 
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Table 4. Evaluation indices of each NLME model. 





Effects 


Mixed parameters 


AIC 


BIC 


-2LL 


Nested effects of plots and trees 


h 


2992.3380 


3022.0250 


2980.3380 




<h 


not converge 










3649.1130 


3678.8010 


3637.1140 




§\< $2 


2164.7940 


2214.2730 


2144.7940 




fy\t <^3 


2080.9480 


2130.4270 


2060.9480 




t l , 2' <^3 


not converge 








^1/ $2/ <j>3 


not converge 






Plots effects 




5151.7190 


5176.4590 


5141.7200 




§2 


5744.9100 


5769.6500 


5734.9100 




h 


5209.6700 


5234.4100 


5199.6700 




^lf <^2 


5146.6480 


5181.2840 


5132.6480 




^1i ^3 


5145.6490 


5180.2850 


5131.6500 




^2/ <^3 


5168.1260 


5202.7620 


5154.1260 




ij>v 4>2, h 


not converge 






Trees effects 


h 


2995.6300 


3020.3690 


2985.6300 






4145.6660 


41 70.4060 


4135.6660 






3651.0750 


3675.8150 


3641.0760 




<l>2 


2167.8120 


2202.4480 


2153.8120 




4>1. <I>1 


2083.0270 


2117.6630 


2069.0280 




<h h 


2102.0420 


2151.5210 


2082.0416 




h, <l>2, h 


not converge 
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AR(1) = CT 2 
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(9) 



ARMA(1,1) = <7 2 



1 

V 

yp y 



y yp 

1 y 
1 



(10) 



cs= 



a 2 + o\ 



(11) 



where p is the autoregressive parameter, y is a moving average 
component, a" is the residual variance, and a l is the residual 
covariance [37-40]. 

Parameter estimation. The parameters in the equations 
were estimated by maximum likelihood (ML) using the Lindstrom 
and Bates (LB) algorithm implemented in the R nlme function 
[18,22]. The LB algorithm and nlme function are detailed in 
several articles (see, for example, [18,22]). 

Predicting the random effects parameters is more problematic 
during model application and prediction than during the fitting 
process. In this case, they were estimated by the EBLUPs [25], 
using the increment data. 



b.xDZt R,+Z,DZM e, 



(12) 



where bi is the estimated random effects vector of EBLUPs, D is 
the qxq estimated variance-covariance matrix (q is number of 
random-effects parameters) for the random effects, R/ is the 
estimated variance-covariance matrix for the error term, Z,- is the 
estimated partial derivatives matrix with respect to the random 
effects parameters for the new observation, and e, is the residual 
vector, whose dimension is the number of observations, and whose 
components are given by the difference between the observed 
diameter growth value for each tree, and the value predicted by 
the model including only fixed effects. 

The standwise calibration was used to evaluate the accuracy of 
the calibration [12]. This type of calibration involves using the 
random plot components predicted from the increments of a small 
sample of trees per plot to predict the increment of the trees within 
the plot not used in the calibration. In this case, the calibration was 
made with 1, 2 and 3 trees per plot. The random parameters of 
new observations could be predicted with Equation 12. 

Results 

Function selection 

The R Tils function was used to evaluate the parameter 
estimates and model fit statistics of the five equations (Table 2); the 
results are listed in Table 3. The Korf equation had slightly better 
predictive ability than the others. Therefore, the Korf equation 
was selected as the basic nonlinear model for estimating diameter 
growth. The final base model is given by: 
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Table 5. Performance criteria for the best NLME models. 







Effects 


Mixed parameters 


AIC 


BIC 


-2LL 


LRT 


p Value 


Plots effects 


fa 


5151.72 


5176.46 


5141.72 








fa, fa 


5145.65 


5180.29 


5131.65 


10.07 


0.0065 




fa, fa with exponential 
function and ARMA(1,1) 


2042.40 


2091.88 


2022.40 


3109.24 


<0.0001 


Trees effects 


fa 


2995.63 


3020.37 


2985.63 








fa, fa 


2083.03 


2117.66 


2069.03 


916.60 


-C0.0001 




fa, fa with exponential 
function and ARMA(1,1) 


1113.96 


1163.44 


1093.96 


975.07 


<0.0001 


Nested effects of 
plots and trees 


fa 


2992.34 


3022.03 


2980.34 








fa, fa 


2080.95 


2130.43 


2060.95 


919.39 


<0.0001 




fa, fa with exponential 
function and ARMA(1,1) 


1112.75 


1177.07 


1086.75 


974.20 


<0.0001 
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DBH 



ijk - 



-<i>i exp 



* //a 



-B ijk 



(13) 



NLME model construction 

The approach used to construct the NLME models was to fit 
the models with nested effects of plot and tree for Equation 13 and 
then to successively remove the random effects. The results are 
listed in Table 4. Four of the NLME models reached convergence 
with nested effects of plot and tree; the fifth and sixth models 
converged when the one-level models included the random effects 
of plot and tree, respectively. 

LRT, AIC, BIC and —2 LL fit statistics were compared among 
different combinations of random effects parameters (Table 4). 
The models represented by Equations 14—16, incorporating the 
nested effects of plot and tree, plot effects and tree effects on <I> 1 
and (j> 3 , yielded the smallest AIC, BIC and — 2 LL. 



DBH? (t and 



=(/f l +«i,+w 1 , y )exp -ft / q k 



+ 8* (14) 



NLME models with heteroscedasticity and 
autocorrelation 

We used the power function or the exponential function as the 
variance functions and the AR(1), ARMA(1,1) or CS as the 
autocorrelation structures to fit diameter growth models incorpo- 
rating different random effects. The results of the models provided 
the best fit are shown in Table 5. The selected models had the 
smallest AIC, BIC and —2 LL. Thus, the final models of plot 
effects, tree effects and the two nested effects are, respectively: 

Equation 14 + Equation 7 + Equation 10 (17) 



Equation 15 + Equation 8 + Equation 10 (18) 



Equation 16 + Equation 7 + Equation 10 (19) 



Parameter estimates 

Nested effects of plot and tree. The nested two-level 
NLME diameter growth model is: 



DBHP, = (/? 1+M| ,)exp( -h/tfc+^+B* (15) 



DBH* -(ft+HyJrap ~h % 



-Bijk (16) 



where 



DBHP, and ' = (27.2796 



exp -24.7556 



+ U\i +Uiij) 



( i.5416+«,,-+h,,,)\ (20a) 
4 



where DBH^ and 1 and DBH? /f , 



DBHj /Ar are the diameters at 



breast height for the three effects; j8 2 and /J 3 are fixed-effects 
parameters; Uu and W3, are random-effects parameters generated 
by plot on 0j and <f> 3 , respectively; Uy and Uy are random-effects 
parameters generated by tree on 0j and </> 3 , respectively; and U\y 
and U3ij are random-effects parameters generated by interaction of 
plot and tree on (f>i and ^3, respectively. 



/'/ = 



~u u ~ 




0~ 






0_ 









Ulij 



,^ij : 



' 3.5713 
v -0.1411 

6.0658 
-0.5052 



-0.1411 
0.0640 



-0.5052 V 
0.0788 J 



(20b) 



(20c) 
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Equation 1 


3 




■4ff| 




''lijjiiil;:::::^ 



Equation 20a 



0 5 10 15 20 25 

Predicted values(cm) 



10 15 20 25 
Predicted values(cm) 




Equation 22a 



10 15 20 25 30 35 
Predicted values(cm) 



0 5 10 15 20 25 30 35 
Predicted values{cm) 



Figure 2. Residual error map of diameter growth of each model. 
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Sijk 



~N(0,R ijk \.iPi)i)<n..(r y ) (20d) 



varexp (s iJk ) = 1 .0700exp ( - 0.0724%) (20e) 



ARMA(1,1) = 1.0700 x 



1 0.7842 0.2456 
0.7842 1 0.7842 
0.2456 0.7842 1 



(20f) 



Plot effects. The NLME diameter growth model incorpo- 
rating the effect of plot is: 



DBHP. t =(21.9361+« 1 Oexpf-52.3880 A^' 9153+ " 3,) )+%/ f (21a) 



where 



5.2409 
-0.2937 



-0.2937 
0.0910 



(21b) 



ARMA(1,1) =2.3484: 



1 0.9822 0.1811 
0.9822 1 0.9822 
0.1811 0.9822 1 



(21e) 



Tree effects. The NLME diameter growth model incorpo- 
rating the effect of tree is: 



DBHk=(27.2695 + «y) 
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Figure 3. Scatter plot of fitted values against observed values of diameter growth of each model. 
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Model prediction 

The predictive ability of Equation 13 was evaluated using 
predict procedures and Equations 3-5 on both fitting and 
validation data. The performance of the NLME models, with 
and without modeling the error structure, was evaluated using 
cross-validation procedures for both fitting and validation data; the 
random effects were predicted with the EBLUPs (Equation 12), 
using the measurement data. 

Table 6 lists the three fit statistics for Equation 13 and 
Equations 20a— 22a with and without random effects. Equation 
20a was the best predictor, with increases in R^ and decreases in 
AMR and RMSE for both fitting and validation data, but it was 
more complex than the others and incurred significant computing 
cost. In Figure 2, the residuals of Equations 13 and 20a-22a are 
plotted against the fitted values; the fitted values are plotted against 
the observed values in Figure 3. Based on the above analysis, we 
can conclude that, although Equation 20a is the strongest 
predictor, it is more complex than Equation 22a and the 
difference between them is small. Compared with Equation 13, 
Equation 22a had a higher R\q, 0.9956 compared to 0.7758, and 
a lower RMSE, 0.5344 compared to 3.2810. Therefore, the 
NLME model incorporating the random effect of trees was the 
best model for predicting diameter growth of individual China-fir 
trees in the single-species plantations of the study area. 

Discussion 

Of the 5 theoretical growth equations tested, the Korf equation 
best fit the individual-tree diameter growth data of China-fir when 
evaluated on the basis of AMR, RMSE and Rl dj . The Korf 
equation is widely used for forest growth and yield simulation 
models [41—43]. ONLS regression is commonly used to build 
forest growth models, but its value is limited because tree data 
typically violate the assumption of independent and identically 
distributed errors [13,44,45]. NLME models are a useful tool for 
analyzing repeated measures data and spatially correlated data 
[18,33]. A model can be constructed with a unique variance- 
covariance structure that eliminates the influence of the random 
effects (plot and tree effects in this study). The two primary 
challenges in fitting NLME models are determining the mixed 
parameters and calculating the random effects [9,18,24,33,46]. An 
additional source of inherent correlation would be the effect of 
year, where observations coming from the same year would be 
highly correlated; tree-ring width is largely related with yearly 
climate variables [47]. However, year effects were not analyzed in 
this study. Incorporating annual climate factors into the NLME 
models may be an appropriate area for future research. 

The Korf equation has been widely used as the base NLME 
model for forest growth and yield prediction. For example, Cheng 
and Gordon [48] successfully used the Korf equation with NLME 
models to fit loblolly pine (Pinus taeda L.) diameter-age 
relationships; the one-level (tree) individual-tree NLME model, 
based on the Korf equation, with random effects parameters (j>i 
and (^3 had the best fit. Parameters (j>i, <j> 2 and ^ 3 are the 
asymptotic values, the values associated with the growth rate of the 
tree and the values associated with the curve shape (inflection 
point) of the Korf equation, respectively. Therefore, the random 
effects (tree) mainly influence the maximum value and the 
inflection point, with evidence that the growth rate of the tree 
affects the model fit. 

Sometimes, no prior information is available from which the 
random parameters can be predicted. In this case, the mixed- 
effects model with the random parameters set to 0 is not the same 
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as the population average model and will give biased predictions. 
Instead, the population average model, fit without random effects, 
should be used. 

Conclusions 

Five theoretical growth equations were evaluated for estimating 
the diameter growth of China-fir trees grown in monospecific 
plantations in Fujian province, southeast China. The equations 
can be evaluated for both biological and statistical meaning. All 5 
equations and the Korf equation in particular were commonly and 
successfully used to model individual-tree diameter growth. One- 
level (plot or tree) and nested two-level (tree nested within plot) 
NLME models based on the Korf equation, with variance 
functions and correlation structures, were used to estimate 
diameter growth of individual trees; this approach was necessitated 
by the hierarchical structure of the experimental design and the 
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